Photo-induced separation of chiral isomers in a classical buffer gas 
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We develop a theory of photo-induced drift of chiral molecules or small particles in classical 
buffer gases. In the absence of a magnetic field there exists a flux of chiral molecules, provided the 
electromagnetic field is circularly polarized. It has opposite signs for different chiral isomers. In the 
presence of a magnetic field the flux can be also induced by a linearly polarized (or unpolarized) 
electromagnetic field. The magnitude of the flux is not proportional to either linear or orbital 
momentum of the electromagnetic field. 
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In this article we develop a theory of photo-induced 
transport of chiral molecules (or small chiral particles) 
diluted in a classical buffer gas. This effect can be impor- 
tant for separation of different isomers of chiral molecules 
and for the theory of optically actuated molecular motors 
(see for example [H G2 H El ) • 

Consider a solution of molecules of definite chiral- 
ity in a buffer gas in the presence of an electromag- 
netic radiation with frequency u and polarization tensor 
Pij = (EiEj), where E is electric field and (...) denotes 
time averaging. We define the chiral current j c as the 
part of molecular current that changes sign when the 
molecules are replaced by stereo-isomers of opposite chi- 
rality. We study effects whose magnitudes are not pro- 
portional to either linear or angular momentum of light, 
i.e. we work in the approximation where no force or 
torque is exerted by the electromagnetic radiation on the 
medium. This approximation is valid as long as the ther- 
mal angular momentum of the molecules is larger than 
h. In this case in the absence of a magnetic field, by 
symmetry, the chiral flux density is nonzero only if the 
electromagnetic field is circularly polarized, 
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Here is the antisymmetric tensor, and the coefficient 
a c has opposite signs for molecules of different chirality. 
It is clear from Eq. that when the direction of circular 
polarization is reversed the chiral current changes sign. 

Although the photo-induced flux of chiral molecules ex- 
ists both in gases and liquids the microscopic mechanisms 
of the effect in these two cases are different. In this arti- 
cle we focus on the case of gases, where all relevant mean 
free paths (for both buffer gas and chiral molecules) are 
bigger than sizes of the molecules. In this case the effect 
can be analyzed using the Boltzmann kinetic equation. 

We assume that optical transitions take place be- 
tween the ground state and an exited state of the chiral 
molecules denoted by indices g and e respectively, the 
energy difference between them being Hujq. 

Let us first illustrate the origin of the effect by con- 
sidering a toy model in which the chiral molecules are 
constrained to rotate about a single axis which coincides 
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FIG. 1: (a) Schematic drawing of the dependence of the dis- 
tribution function of chiral molecules on their angular velocity 
along the wavevector of light. 



with the direction of propagation of light. The frequency 
of the circularly polarized radiation in the reference frame 
co-rotating with the molecule is shifted by the rotation 
frequency of the molecule, Q. As a result the proba- 
bilities of light absorption by molecules rotating clock- 
wise and counterclockwise are different. Consequently 
the distribution function of the molecules becomes asym- 
metric with respect to the angular velocity: a hole will 
occur in the distribution function of the ground state 
molecules f 9 (Q) at certain values of the angular velocity, 
while a peak will occur in the exited state distribution 
f e (Q)(see Fig. [I]). Thus the populations of molecules in 
the excited and ground states acquire opposite angular 
momenta. A fundamental property of collisions of the 
chiral molecules with the buffer gas is that they trans- 
fer rotation of chiral molecules into translation. Since 
the scattering cross-sections for the excited and ground 
states are generally different this will result in the net flux 
of the chiral molecules. Momentum conservation during 
collisions implies that a net momentum in the direction 
opposite to the flux is imparted to the buffer gas. In some 
aspects this effect is similar to the light-induced drift of 
atoms considered in Refs. [6j [7J |9j 

We now turn a description of the effect in the frame- 
work of the Boltzmann kinetic equation. Due to the rapid 
mixing of rotational degrees of freedom of the molecules, 
their distribution function depends only on the center of 



mass coordinate r, time t, and the set of quantities T that 
are conserved during free motion [10] . Below we consider 
a spatially-uniform situation and assume that the chiral 
particles scatter only from non-chiral buffer gas, whose 
distribution function is in equilibrium. Then the Boltz- 
mann kinetic equation has the form 

= /°{/(r)} + /-„(r). (2) 

Here / a (r) is the nonequilibrium part of the distribution 
function, and a = g,e labels the ground or exited states of 
the chiral molecules. Denoting the scattering probability 
of a molecule from state a, T to a' , T' by w a a (r / ;T) we 
write the scattering integrals in the excited and ground 
state due to collisions with the buffer gas as 

I 9 {fF)} = J w^(T;T')fS(r')dr' 
-f 9 (T) J w"(T';T)dr' 
+ J w° e (T;T')f(T')dT', (3a) 

/ e {/(r)} = -f(r) j w^{Y'-Y)dY'. (3b) 

To linear order in the light intensity X the transition rates 
between the ground and exited states due to light absorp- 
tion are described by 

#.(r) = -/o(r)T(r)j, (4a) 
i; h (r) = / (r)T(r)j, (4b) 

where /o(r) is the equilibrium distribution function of 
the molecules in the ground state (we assume that the 
excited state is not thermally populated), and Y(T) is 
the probability of light absorption. The latter is deter- 
mined by the time dependence of the electric field in the 
reference frame of the molecule, and in general depends 
both on the molecule velocity (due to the Doppler effect) 
and its angular velocity. The rotational frequency shift 
is of order of the angular velocity of rotation, ~ vr/d, 
where vt is the thermal speed of the molecule and d is 
its characteristic size. The Doppler frequency shift is 
much smaller, of order vt/\ where and A is the wave- 
length of light. We assume that the mismatch between 
the radiation frequency uo and the resonant frequency 
ujq is of the order of the thermal rotational shift, i.e. 
\uj — uoq\ ~ Vt/cI, and that the absorption line width 7 
is smaller than the rotational broadening 7 < iv T /d. In 
this regime the Doppler effect can be neglected. 

Thus we take the absorption coefficient T(T) to 
be independent of the linear momentum of the chiral 
molecules. Since we work in the dipole approxima- 
tion, in which T(T) is insensitive to the molecule chi- 
rality it is invariant under inversion. In other words 
I£ h (r) = Ip h (r p ), where the set of coordinates T p is 



obtained from T by inversion. If the molecules were non- 
chiral the scattering probabilities would be invariant un- 
der inversion, i.e. 

w(T';T)=w(r' p ;T p ). (5) 

In this case the nonequilibrium part of the distribution 
function f(T) would be a symmetric function of momen- 
tum, /(r) = f(T p ). This would nullify the chiral current 

j c = £ / drv /°( r ), (6) 

where v is the molecule velocity. For chiral molecules, 
in addition to inverting velocity, the inversion symme- 
try changes the molecule chirality. In this situation the 
equality Eq. ([5| is not valid, and the nonequilibrium dis- 
tribution /(r) contains a part that is odd in momentum 
and the chiral current Eq. (|6| generally does not vanish. 

In order to complete the description of kinetics of the 
molecules in the presence of radiation we need to ob- 
tain the dependence of the absorbtion probability on the 
phase space coordinates, Y(r). In the dipole approxi- 
mation T(r) can be expressed in terms of the three fre- 
quency dependent absorbtion cross-section, ai(u), for the 
different light polarizations, by transforming the electric 
field from the lab frame to the reference frame co-rotating 
with the molecule. This problem can be solved for the 
general case, where molecule rotation is described by an 
asymmetric top [TT] . However, in order to simplify the 
discussion and illustrate the essential phisics we consider 
the case where the molecule rotation is modeled by a 
symmetric top. 

In this case the set of conserved quantities consists 
of the linear momentum p, angular momentum M, and 
the angle between the molecular axis and M, i.e. 
T = {p, M, 6} (we introduce the Euler angles 0, 0, and 
ip, as shown in Fig. [2] b))[20]. In the case of circularly 
polarized light the distribution function depends not on 
all three components of the angular momentum M but 
only on its magnitude M and the angle (3 it makes with 
the direction of light propagation (z axis in our nota- 
tions, see Fig. [5] a)). Thus the distribution function / 
depends only on p, M, /3. Furthermore, it is conve- 
nient to express M and in terms of the components 
of the angular momentum parallel, Ms = Mcos#, and 
perpendicular, M± = Msin#, to the symmetry axis of 
the molecule. Indeed the rotational energy of a sym- 
metric top is most naturally expressed in these variables; 
e(T) = 277* Therefore below we write the equi- 

librium distribution function and the kinetic equation in 
the variables T = {p, M3, M_l, /?}, with the integration 
measure dT ~ d 3 p d cos /3 M±dM±dMs. 

We assume that the excited state is non-degenerate. In 
this case, for a symmetric top molecule and in the dipole 
approximation, only the electric field polarized along the 
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FIG. 2: (a) Rotation by the angle (3 from the xyz frame to 
the frame, (b) Illustration of the Euler angles 0, <fi, ip 
rotating the £77^ frame to the x 1X2X3 frame co-rotating with 
the molecule. 



symmetry axis, £3, of the molecule is absorbed. Further- 
more, the absorption cross-section for the right and left 
stereo-isomers are identical and will be denoted by a(cu). 

In order to express T(r) in terms of a(u) we need to 
transform the electromagnetic field E from the stationary 
frame to the frame rotating with the molecule. To this 
end we introduce three reference frames as illustrated in 
Fig. [2j The xyz frame is the lab frame such that the 
z-axis corresponds to the direction of light propagation. 
The light is assumed circularly polarized with the polar- 
ization vector Ei = Eo(x -\- iy) / \/2, where the hat above 
a symbol is used to denote a unit vector along the corre- 
sponding axis. The £77^ frame is also stationary but the 
£ axis points along the constant angular momentum M 
of a given molecule and forms an angle (3 with the z-axis. 
The £ and the x axes coincide. In terms of the basis 
vectors £, 77, £ of this frame the electric field is written as 

with the components expressed in terms of Eq as 
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Next we transform the electric field from the fixed frame 
^r]( to the frame X1X2X3 which rotates with the molecule 
so that the X3-axis always point along the symmetry axis 
of the molecule. The rotation from the £77^ frame to the 
X1X2X3 frame is described by the Euler angles (p,0,ip. 
The components of the electric field in the rotating 



frame, E = £ + *i±p + E x 3 + E 
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are given by 
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For a symmetric top molecule the angles <p and ip change 
linearly with time while the angles (3 and do not change 
during free propagation. Thus it is clear from Eq. ^ that 
in the reference frame co-moving with the molecule each 



component of the electric field, E±, and Eq, consists of 
three monochromatic waves with frequencies of the form 



u pq = u) - pip - q<\>, 



(9) 



where p and q can take the values p,q = 0, ±1. Since 
only the field £"0 along the molecule axis X3 is absorbed, 
we can write the absorption probability using Eq. |8) as 
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The frequencies uo pq can be expressed in terms of the 
variables M 3 and Mj_, M = y^M' 
matic relations 
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Equations ( [To] ) ? (JqJ) and (11) define the absorption 
probability as a function of the phase space coordinates 
T = {p, M 3 , M_l, (3}. This probability is peaked when 
the resonance condition pip + q<p = uo — uoq is satisfied. In 
terms of the variables M3, Mj_ the resonance condition is 
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■Ml =I ± (lj-lj ). 



We note that the light absorbtion probability Eq. (10) 



contains a term that is linear in cos (3. The coresponding 
optical transitions produce a nonequilibrium distribution 
function that is odd in the component of the angular mo- 
mentum along the light propagation direction, M z . Sub- 
sequent collisions of the chiral molecules with the buffer 
gas result in a net flux of the molecules along the z-axis. 

In order to obtain an estimate for the magnitude of the 
effect we consider the simplest case when the relaxation 
time r eg from excited to ground state is the shortest time 
in the problem. Here l/r ge (T) = / dT'w 9e (T, T'). Then 
in a stationary state we get from Eqs. ([2]), ( |3b[ ) and (4b) 



/ e (r) = T» e (r)rr(r)/ ff (r), 



(12) 



while the Boltzmann equation for the ground state dis- 
tribution function becomes 

f 9 (T) [ drv^(r' ; r)- / dr'w 99 (T;T , )f 9 (r , ) = 



-T(r)/#(r) + J dT , w 9e (T',T , )r 9e (T , )T(T , )f 9 (T , ) 



In our approximation the molecules in the excited state 
do not contribute to the total flux of chirality, i.e. 



/ 



vf 3 (T)dT. 
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Thus we are interested only in the anisotropic in v part 
of f 9 (T), which can be estimated as 

j c ~ r 99 i J drdr'vr(r')f°(r') x 

[ w 9z(T,T f )T 9e (T f ) -w 99 (T,T f )r 99 (T% (13) 

Here we introduced the relaxation time for scattering of 
chiral molecules in the ground state, r" , and the mo- 
mentum relaxation time, r 99 . 



A rough estimate for the chiral current Eq. (13) can 
be written as 



j c ~ v T n c (Br p )(5W 9e r 9e - SW"r") 
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(14) 



Here n c is the concentration of the chiral molecules, 
vt is their thermal velocity, £ is the absorption length 
of the light, and B ~ n' 1 ! J dTT(T)fo(T) is the rate 
of the optical transitions in an individual molecule. 
Here we introduced the "chiral part" of the scattering 
rate as SW ga ~ (Bu^t)- 1 f dTdV'z ■ v[w 9a (T;T f ) - 
w 9°(FP ; r' p )]fll(r')r(r')l, with a = e/g. Thus 
SW 9a r ga is a dimensionless measure of the degree of chi- 
rality of the molecules (we are not aware of any studies 
of this quantity). Equation (14) holds at small enough 
intensities of radiation, when / <C fo- At large intensities 
the magnitude of the effect saturates, and the maximum 
flux of chiral molecules is of order 
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(J) (SW ge -5W")t 99 . (15) 



We note that relative signs of SW ge and SW" , as well as 
the sign of the chiral current, are arbitrary and depend 
on the structure of the molecules. 

The effects considered above are direct analogues of 
photo-galvanic effects in non-centrosymmetric crystals. 
There are two types of photo-galvanic effects. The first 
one is associated with a transfer of momentum from light 
to electrons, and its magnitude is proportional to the 
momentum of light. The analogue of this effect, which is 
associated with transfer of angular momentum of electro- 
magnetic field to molecules in classical gases and liquids 
has been considered in Refs. [12j [13j [14j [15j The sec- 
ond type of the photo-galvanic effects exist only in non- 
centrosymmetric crystals and is not proportional to the 
momentum of light p~6j [T71 [18j [19] . Since the momentum 
of light is small the magnitude of this effect is generally 
bigger than that of the first one. The photo-induced drift 
of chiral molecules diluted in a buffer gas considered here, 
is an analogue of the second type of photo-galvanic effect. 

There exists also an analogue of the linear photo- 
galvanic effect [I6j El [THJ [19] . Namely, in the presence of 
a magnetic field a flux of chiral particles exists even if the 
electromagnetic field is linearly polarized or unpolarized, 



The photo-induced chiral drift also exists when chi- 
ral molecules are dissolved in a liquid. Photo-induced 
propulsion of molecules in this situation can be described 
in the framework of hydrodynamics. These two aspects 
of the phenomenon will be considered in a separate pub- 
lication. 
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